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The energy-momentum conservation laws for general reduced-fluid (e.g., gyrofluid) models are 
derived by Noether method from a general reduced variational principle. The reduced canonical 
energy-momentum tensor (which is explicitly asymmetric and has the Minkowski form) exhibits 
polarization and magnetization effects associated with dynamical reduction. In particular, the 
asymmetry in the reduced canonical momentum-stress tensor produces a non-vanishing reduced 
intrinsic torque that can drive spontaneous toroidal rotation in axisymmetric tokamak plasmas. 

I. INTRODUCTION 

Nonlinear reduced-fluid models play an important role in our understanding of the complex dynamical behavior of 
strongly magnetized plasmas. These nonlinear reduced-fluid models, in which fast time scales such as the compressional 
Alfven time scale have been asymptotically removed, include the reduced magnetohydrodynamic equations [H-0 , the 
reduced Braginskii equations [j, [B[ , the nonlinear gyrofluid equations @, 0] , and several truncated reduced- fluid models 
(such as the Hasegawa-Mima equation and the Hasegawa-Wakatani equations [13] ). Because the space-time-scale 
orderings for these reduced-fluid models are compatible with the nonlinear gyrokinetic space-time-scale orderings 0] , 
they provide a very useful complementary set of equations that yield simpler interpretations of low- frequency turbulent 
plasma dynamics in realistic magnetic geometries. 

The role of plasma flows in self-regulating anomalous transport processes in turbulent axisymmetric magnetized 
plasmas has been intensively investigated in the past decade. Because a strong coupling has been observed [III Ell 
between toroidal angular-momentum transport and energy transport in such plasmas, it is natural to investigate the 
link between these two global conservation laws through an application of the Noether method on a suitable Lagrangian 
density. Axixymmetric tokam ak p lasmas have also been observed to undergo spontaneous toroidal rotation in the 
absence of external torque [l3l - ll7l |. Several attempts have been taken to explain this surprising phenomena [l8l - [23| . 
In the present work, we apply the Noether method on a general reduced Lagrangian density to identify an intrinsic 
torque generated by the dynamical reduction, which introduces a decoupling between the reduced-fluid momentum 
and the reduced-fluid velocity as well as reduced polarization and magnetization effects. 

The remainder of this paper is organized as follows. In Sec. lU we present the variational formulation of general 
reduced-fluid equations that exhibit polarization and magnetization effects in the reduced Maxwell equations as well as 
a reduced ponderomotive force in the reduced-fluid momentum equation. In Sec. Mil we present the Noether derivation 
of the reduced energy-momentum conservation laws that are preserved by the reduced-fluid equations. Explicit proofs 
are also presented that identify the energy-momentum fluxes as well as the energy-momentum transfer terms. The 
reduced toroidal angular-momentum conservation law is presented in Sec. IIV1 where the source of reduced intrinsic 
angular-momentum is identified as the reduced intrinsic torque due to the antisymmetry of the reduced momentum- 
stress tensor. We show how this asymmetry is explicitly due to the decoupling of the reduced-fuid momentum and the 
reduced-fluid velocity as well as the polarization and magnetization effects associated with the dynamical reduction. 
We summarize our work in Sec. |V] and discuss future work. 

II. GENERAL VARIATIONAL FORMULATION 



The general variational formulation of nonlinear dissipationless reduced-fluid models describing nonlinear turbulent 
dynamics of strongly magnetized plasmas is expressed in terms of a Lagrangian density C(ip a ) that is a function of 
the multi-component field 

V" = ($ J A ) E,B;n,u,p||,pj.) J (1) 
where the electromagnetic fields (E, B) arc defined in terms of the electromagnetic potentials ($,A) as 



E = — V3> — c^dA/dt and B = V X A 



(2) 
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and the reduced-fluid moments (n,u,pt\,p±) are used for each plasma-particle species (with mass m and charge q). 
Gauge invariance requires that the potentials ($,A) should not appear in the final form of the reduced dynamical 
equations as well as their associated conservation laws. 

We note that the Lagrangian formalism does not accommodate higher-order fluid moments (e.g., heat fluxes) and, 
therefore, the issue of fluid closures is completely ignored (unless these higher-order moments are expressed in terms 
of lower-order moments). These higher-order moments, as well as dissipative effects, can thus be added after the 
dissipationless reduced-fluid equations are derived by variational method (although a variational procedure (2~3 . |25| 
can be used to include heat fluxes in the pressure evolution equations). The emphasis of the present work is therefore 
on nonlinear physics within a Lagrangian formulation of reduced-fluid dynamics. 



A. Reduced-fluid Lagrangian 



We now introduce the general Lagrangian density 

C{^p a ) = £ M (E,B) + £*($,A;n,u) + £ F (n,u,p||,p ± ;E,B), (3) 
defined as the sum of the electromagnetic Lagrangian density 

£ M (E,B) = ^(|E| 2 - |B| 2 ), (4) 

the gauge-dependent plasma-electromagnetic interaction Lagrangian density 

£*($,A;n,u) = - ^ - A- = - ^ qn*, (5) 

where ^2 denotes a sum over particle species, and the reduced-fluid Lagrangian density Cf, which depends on the 
electromagnetic field (E, B) only through the process of dynamical reduction [jt| . 

Next, we introduce the following partial derivatives of the Lagrangian density First, we introduce the 

electromagnetic-potential derivatives 

which define the reduced charge density g = ^ qn and the reduced current density J = ^ qn u. The definitions © 
ensure electromagnetic-gauge invariance, which is connected to the charge conservation law 

d (dC\ „ ( dC\ dg „ _ 

as follows. Under a gauge transformation $' = $ — cr^dxjdt and A' = A + V%, we find ty' = $ — c^dx/dt, where 
d/dt = d/dt + u • V , so that the Lagrangian © transforms as 

^ £ + 4 (f) + J . v(f)s£ + ! H) + v . (jf) , 

where Eq. (J7]) was used. Hence, the reduced variational principle 

C(^ a ) d 4 x = (8) 



is gauge-invariant since J C d A x = J C d A x. We therefore see that gauge invariance and charge conservation are 
intimately connected within a Lagrangian formulation. 

Second, we introduce the electromagnetic-field derivatives [26M28| 



dC dC\ /D H\ /E „ B . 

9E' 9B ) S W-to) S U +P '^ +M) - W 



which define the reduced polarization and magnetization 

V 8F, ' dB 
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These definitions correspond to polarization and magnetization effects associated with reduced electric and magnetic 
dipole moments [26, 28]. Higher-order multipole contributions require the reduced-fluid Lagrangian density to depend 
on gradients of the electromagnetic fields (which are not considered in the present paper). We note that gyrofluid 
models @ and gyrokinetic models Q implicitly contain all multipole polarization and magnetization contributions 
whenever finite-Larmor-radius effects are retained to all orders. 
Third, we introduce the reduced-fluid derivatives 

K = dC F /dn ) 

\, (11) 
p = n 1 dCp/du J 

which define the reduced kinetic energy K and the reduced kinetic momentum p. Lastly, we introduce the symmetric 
reduced pressure tensor 

P* = -2 P ||^bb - ( r -bb) = P||7||bb + p± 7 ± (i-bb). (12) 

We note that standard finite-Larmor-radius (FLR) corrections are automatically included in the reduced pressure 
tensor flT2J| through the term 7j_ ^ 1 [29j . 

reduced-fluid models derived from the variational principle © will be distinguished by the physical effects contained 
in the fluid Lagrangian density ([3]) through the fields defined in Eqs. (|6l)- (fT2|) . 



B. Constraint Equations 

The variational principle © used to derive the nonlinear dissipationless dynamical equations relies on certain 
constraint equations that must be preserved by the field variations 5ij> a . First, the electromagnetic variations 

<5E ee - V<5$ - c- x 8SA/dt and SB = V X 5A (13) 

preserve the electromagnetic constraints 

VxE = -c^dB/dt and V • B = 0. (14) 

For the reduced-fluid moments (n,u,p\\,p±) associated with each particle species, the constraint equations are the 
continuity equation (associated with the conservation of particle number) 

= - V- (nu) , (15) 



Of 



and the pressure equations 



V- [p\i uj - 2p|| bb : Vu, (16) 
V- (p ± u) - p ± fl- bb) : Vu, (17) 



dp\\ 
dt 
dp± 
dt 

where the Chew-Goldberger-Low (CGL) pressure tensor is [13] 

P = p|| d b + p ± (l - K b) . (18) 
Recall that the CGL pressure equations (fl6| -([T7 |) arise from the adiabatic conservation laws [30l. [3ll] 

l(ls)=». ^ 
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associated with the adiabatic invariance of the gyro-action (or magnetic moment) and the bounce-action, respectively, 
where dB/dt = b • V X (u X B) is used in Eqs. (fl"9]) - (f20"]) . The variations (Sri, Su, Sp±, 5p») that preserve the constraint 
equations ([TSJ and (fj"9" ]) -([2"0" ]l are 

-V- (n€) 



5n 
Sp± 



(d/dt + u-V)£ - |-Vu 

- V- (p ± €) - P±(I-bb):V€ 



(21; 



(Spy = - V- (p|| €) - 2py bb : V£ 

where £ generates a virtual spatial displacement for a fluid element of each particle species and SB = b • V X (£ X B) 
is used in order to satisfy the CGL adiabatic conservation laws (|T9")) - (j2T)|) . Note here that magnetic field B = B b 
considered in Eq. (1211) involves the general magnetic field and not just the background time-independent magnetic 
field used in our previous work [2t1 |. 



C. Euler-Poincare Equations 



Using the variations (fTB")) and (f2"Tj). the variation of the Lagrangian density C(ip a ) for a general reduced-fluid model 
is expressed in terms of the variation fields (6$>,SA, £) as 



sc = 



8A 



dt 

- E € 



v-r + 5$ 



<9£ 

a¥ 



9E 



a 5£ „ / ac 

■777 ~r — + V- U — 



at 9u 



9u 



az: 



vu-— - [r, a v— 



9u 



c at aE 

5£ 



a?/ a 



Vx aB 



V-P* 



(22) 



where 77" = (n,pu,p±) combines the density and pressure fluid moments (summation over repeated indices is implied), 
and the Noether density A and flux T are 



A 

r 



E* 
E 



dC 

du 



dC 



dC 



du 



dC 



- (5$ — 
5E 



X a~B ' 



Note that the Noether terms dA/dt + V • T in Eq. (|22|) do not contribute to the variational principle { 
appear as a space-time divergence. 

For arbitrary variations ((5$, <5A, £), the variational principle jSJ yields the Euler-Poincare equations 



dC 

a$ 

5£ 



V- 



3E 5 



1 a /ac 



aA 



c at v aE 



+ ™ + v x — , 



dC 



= 

at lau 



V- I u — 

ou 



aB' 
au 



a r, dC 

ary a 



V-P, 



The first two equations yield the reduced Maxwell equations 



V • D = 4tt q, 

_ TT 1 3D 4vr T 

VxH = _j ; 

c ot c 



(23) 
(24) 
since they 

(25) 
(26) 
(27) 

(28) 
(29) 



where we used the definitions ([H])-©. These reduced Maxwell equations can also be expressed in terms of the 
electromagnetic fields (E, B) as 



V-E 



1 aE 4tt 
VXB - - — = — 

c at c 



4tt (q - V-P 
J 



— h c V X M 

at 



(30) 
(31) 
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where g po \ = — V-P denotes the polarization density, J po i = dP/dt denotes the polarization current, and J mag 
cV X M denotes the magnetization current. 

Equation (|27|) yields the time evolution equation for the fluid velocity u, subject to the constraint equations (fl5 
([T7]) . Using the definitions (TTT])-([n]), Eq. (gTJ) becomes 
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nx- +U-V 



E + - x B 

c 



n(Vif - Vu-p) 



P± V71. + y V7|| 



V-P, 



(32) 



where the first term on the right side is the Lorentz force on the charged reduced-fluid, the terms involving WK and 
Vu • p involve the reduced kinetic energy and the reduced kinetic momentum , while the remaining terms involve 
generalizations of the CGL pressure-tensor force (with 7_l, 711 ^ 1). 

For a regular particle-fluid model, where the particle-fluid Lagrangian is £ p f = mn |u| 2 /2 — V, where V = \ Tr(P) = 
P± +Pll/2 and reduced polarization-magnetization effects are absent, the particle-fluid kinetic momentum is p = mu 
and the particle-fluid kinetic energy is K = m\u\ 2 /2, so that Vif = Vu-p. Next, the particle-fluid pressure 
derivatives are 711 = 1 = j± so that the reduced pressure tensor P* = P is simply given by the CGL pressure tensor 
(fT8|) . Equation ([52)) therefore becomes the particle- fluid momentum equation 



dp 

~dt 



E 



u 



x B - n~ L V-P, 



where the CGL pressure-tensor force 

V • P = Vp±_ + pa b • Vb 
is expressed in terms of the pressure anisotropy pa = P\\ ~ P±- 



V- 



III. REDUCED ENERGY-MOMENTUM CONSERVATION LAWS 

The purpose of the present Section is to derive the reduced energy-momentum conservation laws by applying 
the Noether method on the reduced Lagrangian density (J3J. We also present explicit proofs of energy-momentum 
conservation for these reduced-fluid models in order to uncover the energy-momentum transfer processes. The reduced 
toroidal angular momentum conservation law will be considered in the next Section. 



A. Noether method 

After substituting the Euler-Poincare equations (|2~51) - (|2"T1) into the variation (f2"2"j) . we obtain the Noether equation 

v-r. (33) 



The energy-momentum conservation laws are constructed from the Noether equation (|33p based on symmetries of the 
Lagrangian density (0) with respect to space-time translations generated by (cSt, Sx). 
First, we introduce the variations induced by the space-time translations 



£ = Sx — u 6t 

SC = -St (dC/dt - d'C/dt) - Sx ■ (V£ - V£) 

5$ = Sx-E - c^dSx/dt 

SA = Sx x B + c St E + VSx 
where the gauge-dependent term 

Sx = cSt $ — Sx- A 

will be removed below and 

' dt J V 9ip a ' dt dt dip 01 



(34) 



(35) 
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denote space-time derivatives with the multi-component field tp a held fixed. By substituting the variations (|34[) . we 
obtain the Noether components 



A = 



E 



<5x — u St 



dC 
du 



1 

c 

u M ) ■ ¥ 



<5x-E 



du 

1 dSx\ dC_ 
c~dt) 5E 



Sx X B + c St E 4 
, 5£ 



Sx x B 



cSt E 



I 



<9E' 



V5 X ) X 



Sx — u St 
DC 



dB' 



(36) 



(37) 



The gauge-dependent terms in Eqs. (|3"6")) - (|3"7]) can now be removed as follows. First, we note that the right side of 
the Noether equation (|33)) is invariant under the transformation [32|, [33[ 



A 

r 



A 

r 



V-q 

dqjdt 



(38) 



c V X m 



where q and m are arbitrary vector fields. For example, by substituting the following gauge-dependent vector fields 
inEq. (|38]), 



Sx dC , Sy dC 

q = — ™ and m = — — , 
c aE c oB 



we obtain the new Noether components 



A 

r 



E 
E 



Sx 



u <5i ) • np 



n u p 



c 

<9£ 



&xB + cStF, 



D 

47T 



DE 



BH (B H) I 



I 

4tt 



<5x - 

V4-7T 



uSt 



£$r St, 

C^j Sx 



E X H 



(39) 



(40) 



(41) 



In Eqs. (|40|) - ([4"Tj) . the gauge-dependent interaction Lagrangian density = — g$ + A- J/c will be cancelled by the 
specific form ([3|) of the Lagrangian density C through the term SC^ = — St dC^ /dt — 8x- on the left side of the 
Noether equation (|33l) . 

Now that all gauge dependence has been explicitly removed from our Lagrangian formulation of reduced-fluid 
dynamics, we redefine the Lagrangian density © as C = Cm + £f as the sum of the Maxwell Lagrangian density and 
the reduced-fluid density and the Noether equation (|3"3"]l becomes 



v-r, 



(42) 



where the right side is still subject to an additional transformation 

A standard perturbative approach to investigating turbulent dynamics in magnetized plasmas is to separate the time- 
independent nonuniform magnetic field Bo = V X Ao from the time-dependent dynamical plasma-electrodynamic 
fields. Based on this separation, we can substitute 



d'C 



dB dC 



dt 9B 







in the reduced energy conservation law, and the dynamical plasma-electrodynamic fields are able to exchange energy 
among themselves while conserving total (global) energy. 

The nonuniformity of the background magnetic field B implies that the term V'£ ^ in the momentum conserva- 
tion law drives the exchange of momentum (including angular momentum) between the background magnetic field and 
the dynamical plasma-electrodynamic fields [34j |. The Noether Theorem, however, tells us that the total momentum 
is conserved in the direction of spatial symmetry of the background magnetic field. Hence, the total toroidal angular 
momentum is conserved in an axisymmetric magnetized plasma when the intrinsic angular momentum is taken into 
account (see next Section). 
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B. Reduced energy conservation law 



In a time-independent medium, the reduced energy conservation law 

is expressed in terms of the reduced energy density 

£ = np-u + EP + — f|E| r 
and the reduced energy- density flux 



8tt 



IBI 



E 



u I np • u 



dC F 

drf 



P* -u 



4tt 



E X H. 



(43) 



(44) 



(45) 



The energy conservation law (|43[) is normally used to benchmark the accuracy of reduced-fluid numerical codes (see 
Ref. [271 ] for example). For a regular particle-fluid, the standard expressions for the energy density 



ran 



and the energy-density flux 



u. — |u|- 



V 



V 



8tt 



(|E| 



IBI 



P-u 



47T 



E X B 



are directly obtained from Eqs. (144]) - (1451) . 

To prove the reduced energy conservation law (|4"3"|) , we begin with the partial time derivative of the energy density 



d£ 
~dt 



d{np) ^ 



E dT> 

4^ ' ~dt 



B dB 

4^ ' ~dt 



EdC F du. dC F 9E 



dt V 47r dt 47r dt J \^ du dt 9E dt dt 

Upon using Eqs. (fT5")) and the kinetic-energy term on the right side of Eq. (|4l)]) becomes 

,d£ F 



y »M. U = j.e y 

^ dt ^ 



V • Tiup'U - rf u • V 



Next, using Eqs. (fT4")l and the electromagnetic-energy terms become 



E 3D 

4V~d7 



B dB 



= - V 



(^ EXH ) 



a£ F 9B 



4tt at V4tt / sb at 

while using Eqs. P)-([I2]) and l[I5 |) -([rf ) . the last term becomes (using d'X/dt = 0) 

dC F 



(V-Ph.) -u 



JE, 



^ aCp du , 5£ F <9E 9£ f 



<9E at 



at 



E 



V • (7 ? Q u) 



drj a 



P» : Vu 



dC F dB 



dB dt 



By combining Eqs. (j47|) - ((4"9|) into Eq. (|46)) . we obtain the reduced energy conservation law (|43|) . 



(46) 



(47) 



(48) 



«T- (49) 



C. Reduced momentum conservation law 



In a nonuniform medium, the reduced momentum conservation law [3£ 

an can 



at 



V-T = V'£, 



is expressed in terms of the reduced canonical momentum density 



IT 



E n p 



D xB 

47T C ' 



(50) 



(51) 
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which has the Minkowski form [36| . and the reduced canonical momentum-stress tensor 

1 



T = E P * 



Ea «9£f 
77 dn a 



8tt 



(|E| 



IBI 



B • M 



n u p 



1 

4tt 



D E + B H 



(52) 



While the first two terms in the reduced canonical momentum-stress tensor (|52[) are obviously symmetric, the remain- 
ing terms are not in a reduced-fluid model in which (p, D, H) ^ (mu, E, B). We note that, for a regular particle-fluid, 
the standard expressions for the canonical momentum density 



T-rcan 

P 



E 



and the symmetric canonical momentum-stress tensor 



1 

47T 



ran u 



El 



E X B 

47T C 



IBI 



(EE + BB) 



are directly obtained from Eqs. (|51|) -(|52 j) . It is therefore clear that the asymmetry of the reduced canonical momentum- 
stress tensor (|52p is due to the dynamical reduction that led to the reduced-fluid Lagrangian density C F . 

We now prove the reduced momentum conservation law (j50[) by calculating explicitly the partial time derivative of 
the reduced canonical momentum density ([51]) 



dt 



E 



d(np) 
dt 



d /DxB 



dt \ Aire 



First, from Eq. (|3"2"|). we find 

d(np) 
dt 

Next, from Eqs. flU]) and we find 



V • nup + qn ( E 



x B 



rf V 



dC F 
dr/ a 



Vu- 



dC F 
du 



(53) 



(54) 



d /DxB 



dt \ 4tt< 



BH DE 



47T 47T 

E "?(E + ;XB 



BH 

47T 



1 

4tt 



VE D VB • H 



(55) 



Lastly, we add Eqs. (|54 |) - ([55|) with the identity 

dC F 







E Vu 



du 



Vri a 



dC F 
drj a 



VE- 



dC 
d~E 



VB- 



dC 
dB 



(56) 



V'£ - V 



sfO 



El 



|B| 



E ( Vu -p + d 0) + h ( VE ' D - VB ' H ) ■ 



to obtain the reduced momentum conservation law (150 



IV. REDUCED ANGULAR-MOMENTUM CONSERVATION LAW 



The symmetry properties of the reduced canonical momentum-stress tensor (fB"2"]) is physically connected to the con- 
servation of angular momentum [37l l38l ] . i.e., the exact conservation of the total angular momentum requires a symmet- 
ric physical momentum-stress tensor. Because of the obvious asymmetry of the reduced canonical momentum-stress 
tensor (|52[) , the angular- momentum conservation law is discussed separately from the energy- momentum conservation 
laws discussed in the previous Section. 

According to the Noether Theorem, the component of the reduced-fluid momentum density in the direction of a 
spatial symmetry of the unperturbed magnetic field is conserved. In axisymmetric tokamak geometry, for example, 
where the background magnetic field is 



B = Vlp X + B 0v (iP) Vip 



(57) 
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and the background scalar fields are independent of the toroidal angle tp, we obtain 

We proceed with the derivation of the toroidal angular-momentum conservation law by following two different 
approaches. Both approaches identify the antisymmetric part of the reduced canonical momentum-stress tensor 
Ta = (T — T T )/2 as the source of intrinsic angular momentum that needs to be accounted for in expressing the 
exact angular-momentum conservation law. Here, the components of the antisymmetric part of the reduced canonical 
momentum-stress tensor (|52|. 

(T A ) tj = ^e ljk r\ (59) 



1 

2 

are expressed in terms of the reduced intrinsic torque 



nuXp + ExP+BxM, (60) 



where we used the relations © to obtain E X D/47T = E X P and - BxH/47r = BxM. Hence, we see that the 
canonical momentum-stress tensor (|52"|) is not symmetric as a result of reduced polarization and magnetization effects 
and the decoupling of the reduced- fluid momentum p/rau, defined in Eq. (|lip . from the reduced-fluid velocity u. 
The reduced intrinsic torque (|60|) will now be shown to act as a source of canonical angular momentum. 

A. Direct approach 

By using the axisymmetry condition (JS5J), the dot product of Eq. (|5D|) with dx/dip yields the reduced toroidal- 
momentum transport equation 

dU™ n dx /„ T \ „ /„ dx\ ^ „9x 



(V-T)--V. T.g Wx^vg, 



dt dip V / \ dip J dip' 

where T T denotes the transpose of T and the canonical toroidal angular-momentum density is 

dx 

dip 

Using the definition (|5 1 [) for the toroidal canonical angular momentum, we obtain 



n^ an = n-— =l-xxU. (62) 



v-^ D X B dx 

nr-E^ + ^r'v (63) 

where p v = p • dx/ dip denotes the (covariant) toroidal component of the reduced-fluid momentum. Since the mixed 

second-rank tensor \Jdx/dip = Rip — (pR in Eq. (|6f [) is antisymmetric, where R = Vi? = — dip /dip is perpendicular 
to (p = z X R, we rewrite the reduced toroidal angular-momentum transport equation (|6ip as 

9n= an / dx\ dx ~\ 

1 V • ( T • — — I = — T a : V — — = — X -RJ • t = r z , (64) 



dt \ dip ) dip 

where we used the definition (15^)) and t z denotes the vertical-component of the reduced intrinsic torque (1SU1) . This 
equation shows that r 2 acts as a source of canonical toroidal angular-momentum, which can perhaps explain the 
spontaneous rotation of axisymmetric tokamak plasmas [l3l - [l5l | in the absence of external torque. 

B. Noether approach 

We begin our Noether approach to deriving the toroidal angular-momentum conservation law by noting that the 
left side of Eq. (|50l) is invariant under the transformation 



n = n can + v-s 

T s = T - dS/dt 



(65) 
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where the antisymmetric second-rank tensor S is chosen to cancel the antisymmetric part of the canonical momentum- 
stress tensor (|5^|) . i.e., Ta = 05/ dt. Hence, by using the definition for the reduced intrinsic torque, we obtain 
the definition 

Sij = - Eijk <y k (66) 

for the reduced intrinsic angular momentum <r, which satisfies the evolution equation 

— ET = ^nuXp + ExP+BxM, (67) 

which represents a classical version of zitterbewegung [30l . l40l | in which dynamical reduction (e.g., the averaging of 
highly oscillatory motion) introduces intrinsic (spin) angular momentum. 

According to the transformation (|65l) and the definition of the reduced intrinsic angular momentum (166[) . the total 
reduced momentum density becomes 

n = n can + v • s = n can - ]- v x a, (68) 



which satisfies the invariance property [331 ] V • II = V • n can and evolves according to the reduced momentum transport 
equation 

^ + V-T S = V'Z. (69) 

The toroidal the reduced toroidal angular-momentum conservation now becomes 

dn y /_ <9x 



Ts • «- =0, (70) 



dt \ dip 

where the total reduced toroidal angular-momentum density is 

where we used S : Vdx/dip = — a z . By substituting Eq. ([TTj) into Eq. (j6"9")l . we recover Eq. with r 2 = da z /dt 



n v = n • — = n^ an - a z + v • ( s • — ) , (71) 



C. Toroidal angular-momentum conservation in axisymmetric tokamak geometry 

Another useful form of the reduced toroidal-momentum transport equation (|64p [or Eq. (|70p] is expressed as 



a* v ^ "V at 

where (• • • ) = V -1 §{■■■) J d-d dip denotes the magnetic-surface average (labeled by the magnetic flux ip), with 
J = (VV> X Vi? • V<^) _1 denoting the Jacobian associated with the magnetic coordinates (ip, (p) and V = ^ J d-d dip. 
In Eq. (|72[) . the surface-averaged toroidal angular-momentum flux is 

T% = Vi>-T-^ =J2 nuf Pr - ^ {p* E <? + B 1 H *) > ( 73 ) 

where -B^ = B • Vf/> denotes the V'-component of the perpendicular component of the perturbed magnetic field (since 
B o -VV> = 0). 

In the electrostatic limit (B = Bo, E = — V$) and making use of the quasi-neutrality condition ~ 4ir P^, the 
canonical toroidal angular momentum (|63)) becomes 

n r = E + \ 

where we used the identity Bo X dx/dip = V-0, and the toroidal angular- momentum flux (|73l) becomes 
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Hence, in the electrostatic limit, the surface-averaged toroidal angular-momentum transport equation ([721) becomes 



d_ 

m 



1 d_ 



d<S> 

dip 



(r z ), 



(74) 



where the surface-averaged intrinsic torque involves the reduced intrinsic torque component 

t z = z-^nuXp+(ExP+B xM) 
= mnu\\ z • 



u , - U, 



X b 



(75) 



5J makes use of the gyrofluid relations [27] p = mun bo, E X P = (in the 
{mnu\\) bo X Uj_, where U_l denotes the perturbed E x B velocity (i.e., 



The last expression for r z in Eq. 
electrostatic limit), and Bo X M 

the lowest-order perturbed perpendicular reduced-fluid velocity). Equation ([74"]) was recently derived by Scott and 
Smirnov [4l| , without the reduced intrinsic torque, by considering the time evolution of the gyrocenter moment of the 
(guiding-center) canonical toroidal angular momentum mvu bo — (q/c)ip. 

Lastly, we note from Eq. (|75[) that, since the perpendicular reduced velocities uj_ and U_l are typically driven 
by perpendicular gradients of plasma parameters [27| (e.g., electrostatic potential and plasma pressure), the vertical 
component of the reduced intrinsic torque (|75[) appears to predominantly involve the vertical components of these 
perpendicular plasma-parameter gradients. Hence, we expect that the surface-averaged reduced intrinsic torque {t z ) 
might be greatly determined by the up-down symmetry of the axisymmetric magnetized plasma [23} . 



V. SUMMARY 



In the present paper, we have used a variational principle based on the general reduced Lagrangian density §5§ 
to derive general reduced- fluid dynamical equations ([2"5j) - ([2"7[) subject to the constraint equations (IT4"|) - (ITT1) . These 
reduced equations satisfy the energy-momentum conservation laws (|43[) and l|50p as well as the toroidal angular- 
momentum conservation law (|70j) . 

The case of the toroidal angular-momentum conservation law (|70|) is particularly interesting because it identifies 
a reduced intrinsic torque t generated by the dynamical reduction, which generates polarization and magnetization 
effects as well as a decoupling between the reduced-fluid momentum and the reduced-fluid velocity. The vertical 
component r z = da z /dt of the reduced intrinsic torque was shown in Sec. IIVI to be the source of intrinsic toroidal 
angular- momentum a z , which can generate spontaneous toroidal rotation in the absence of external torque. 

Future work will focus on the physical interpretation of the reduced intrinsic torque ([601) , and the reduced intrinsic 
angular momentum <x it generates, as well as investigating the surface-averaged intrinsic torque (t z ) for several reduced- 
fluid models [42T |. Ongoing work (with N. Tronko) is also investigating the reduced toroidal angular- momentum 
conservation law for the nonlinear gyrokinetic Vlasov-Maxwell equations [7|. 
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